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ABSTRACT 

The aim of the present letter is to explain the 'critical behaviour' observed 
in numerical studies of spherically symmetric gravitational collaps of a 
perfect fluid. A simple expression results for the critical index 7 of the 
black hole mass considered as an order parameter. 7 turns out to vary 
strongly with the parameter k of the assumed equation of state p = kp. 

Numerical simulations of the spherically symmetric gravitational collapse of 
a massless scalar field |I] have revealed a kind of critical behaviour with 
the mass of an eventually formed black hole playing the role of an order 
parameter. More precisely, Choptuik found a relation of the form 

MbhHp-pT- (i) 

where Mbh is the mass of the formed black hole, p* is the critical value of a 
parameter p characterizing the 'strength' of the initial configuration and 7 is 
a 'critical exponent', whose value was found numerically to be 7 ~ 0.37. The 
critical value p* is distinguished by the fact that for p < p* the space-time 
stays everywhere regular, whereas for p > p* an apparent horizon is formed 
signalling the formation of a black hole. Subsequently the same type of 
critical behaviour was observed in the collapse of an ideal liquid ('Radiation 
Fluid') || and the collapse of axisymmetric gravitational wave packets ||]. 
In both these cases the critical index 7 turned out to be close to the value 



7 ~ 0.37 found for the scalar field collapse, suggesting a kind of universality 
as known from critical behaviour in statistical mechanics. 

Obviously in an attempt to explain this behaviour a special significance 
should be attributed to the critical solution obtained for p = p*. For the 
case of the 'Radiation Fluid' Evans and Coleman found that this solution 
is asymptotically self-similar, i.e. invariant under a simultaneous rescaling of 
the radial coordinate r and the time t, whereas the critical solution in the 
scalar field collapse exhibits only a discrete self-similarity JT|. Spherically 
symmetric continuously self-similar solutions are comparatively easy to an- 
alyze, because the field equations can be reduced to a system of ordinary 
differential equations. No such simplification is obtained for solutions ex- 
hibiting only a discrete self-similarity. In fact, Evans and Coleman where 
able to determine the asymptotic form of the critical solution in the case 
of the 'Radiation Fluid' integrating the self-similar equations with certain 
boundary conditions to be discussed subsequently. These authors also pro- 
posed to determine the critical index 7 through a linear stability analysis of 
the critical solution. This is precisely the method we shall employ in the 
following. Since only in the case of the fluid model we are able to determine 
the asymptotic form of the critical solution, we shall restrict the discussion 
to this case. 

The considered Eulerian fluid is described by the energy momentum ten- 
sor = (p + p)U^U u — pg^y with energy density p, pressure p and 4- velocity 
Up, adopting the simple equation of state p = kp. The 'Radiation Fluid' 
considered in |2| corresponds to the special case k = 1/3. The field equations 
were already given in 0, however we shall rewrite them in a form suitable 
for the discussion of the self-similar solutions. The spherically symmetric line 
element is parametrized as 



and the 4- velocity is written as U{ = yl — v 2 , U? = v/yl — v 2 in terms of 
the radial velocity v. Instead of p we use the dimensionless combination p = 
^/Anrp. In terms of the logarihmic variables r = ln(— t) and a = lnr— ln(— t) 
the field equations become the autonomous sytem of PDE's (with a prime 
denoting d/da and a dot d/dr) 

1 . 2/1 .0 1 + kv 2 
-A 2 1 - A 2 - 2p 

2 v H 1 - v 2 



A' /A = ^A~ 2 (l -A 2 - 2p^^) (3) 
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B'/B = A~\2A 2 - 1 + (1 - k)p) 



(4) 
(5) 




— (B + v)p'/p + (1 + Bv)v'/{1 - v 2 ) = 
k 

1 -I- Ik 

(B + v)A'JA + vB'/B - ^—^v + 

1 + k 

Bij^p/p + vv/il-v^-A/A) 



(6) 



— (A; + (1 + k)Bv + v 2 )p/p +(B + 2v + Bv 2 )v'/{1 - w 2 ) = 



l + 2Bv + v 2 )A'/A + (1 + v 2 )B'/B - 




+B(v'p/p + (1 + w 2 )w/(l - w 2 ) - 2vA/A) 



(7) 



The equations for self-similar solutions are obtained by dropping the terms 
with r-derivatives. In this case the first and third equation can be combined 
to yield the algebraic constraint 



The general structure of the set of self-similar solutions has been discussed 
in some detail by Bogoyavlenskii || and Ori and Piran||. Solutions with a 
regular origin tend to the fixed point v = p = B = 0, A — 1 for a ^ — oo 
and constitute a one-parameter family parametrized by the central value po 
of the density p. The generic solution with a such a regular origin turns out 
to develop a singularity (shock front) at the 'sonic point' for some finite value 
of a . Only for a discrete set of po-values the formation of the shock front is 
avoided and the solution stays regular for all values of a. Numerical analysis 
shows that these regular solutions can be labelled by the number of zeros 
of the radial velocity field v. The solution described by Evans and Coleman 
is the one with one zero of v. Varying k > one finds that this type of 
solution exists only up to a certain maximal value k max m 0.888, where the 
singular sonic point changes its character. As discussed in [|], |] there are 
three different cases for the nature of this singular point. It turns out that 
two of its eigenvalues degenerate for k = fc max and the singular point changes 
from an attractive node to an attractive focus. 



A 2 = l 



2p(l + (1 + k)v I B + kv 2 
l-v 2 



(8) 
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One may ask why it is this regular self-similar solution that is of relevance 
for the actual collapse. The reason may be that the latter seems to proceed 
without the formation of a shock-front. It is an interesting question, if this 
behaviour is a consequence of the chosen class of initial data or a more general 
feature due to the chosen equation of state. 

The characteristic feature of solutions with initial data close to those of 
the critical solution is, that they first approach the latter, but eventually 
run away from it. In order to describe this run-away we may employ a 
linear stability analysis as proposed by Evans and Coleman 0. In the linear 
approximation we expect an unstable mode ip of the critical solution with a 
time-dependence ip ~ ip>(a,p)e~ UJT with some characteristic frequency u> > 0. 
The requirement of regularity of ip at the origin and the sonic point acts 
as a boundary condition leading to a discrete frequency spectrum. Since ip 
vanishes for p = p*, we may assume ip(a, p) = c(a)(p — p*). Of particular 
relevance is the metric function A, which vanishes at the position of an 
apparent horizon. For p close to p* it takes the form 

A{a,r,p) = A ss (a) - c{a){p - p*)e™ , (9) 

where A ss denotes the self-similar back-ground. Although the back-ground 
solution turns out to have no apparent horizon (i.e. A ss (a) > 0), the expo- 
nentially growing perturbation will eventually lead to a zero of A (assuming 
p > p* and c > 0). Here we have tacitly taken for granted that the lin- 
ear approximation is still valid for perturbations of the same order as the 
background. Let us hence assume A(a h , T h,p) = for some values r h , cr h 
determined by A ss and c. Using = e r7h+Th we obtain 

A 9s (cr h ) = c(a h )(p-p*)— (10) 

leading to the desired relation 

M BH = r h ~ (p-p*)» • (11) 

For the critical index we read off the expression 7 = Xjio. 

Tab. |l] gives the numerically determined values of 7 for various values of 
k. The value 7 = 0.3558 for k = 1/3 compares very well with value 7 ~ 0.36 
given by Evans and Coleman 0. The value k = 0.88 is the maximal one 
for which a regular self-similar solution could be found. The limit k — > is 
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k 


7 


k 


7 


0.01 


0.1143 


0.45 


0.4400 


0.05 


0.1478 


0.50 


0.4774 


0.10 


0.1875 


0.55 


0.5159 


0.15 


0.2251 


0.60 


0.5556 


0.20 


0.2614 


0.65 


0.5967 


0.25 


0.2970 


0.70 


0.6392 


0.30 


0.3322 


0.75 


0.6834 


1/3 


0.3558 


0.80 


0.7294 


0.35 


0.3676 


0.85 


0.7775 


0.40 


0.4035 


0.888 


0.8157 



Table 1: Critical indices 7 for various values of k. 

singular, since the sonic point merges with the origin. Nevertheless 7 seems 
to approach a non- vanishing limit. As can be seen the critical index depends 
rather strongly on the value of k, hence it is not universal. It would be clearly 
desirable to check the predictions made for 7 for values of k different from 
1/3 by actual collapse calculations. 

Acknowledgement! am indebted to P. Breitenlohner, B. Briigmann, 
F. Miiller-Hoissen and A. Rendall for stimulating discussions. 

Note added: While this letter was in preparation a paper || by T. Koike, 
T. Hara and S. Adachi on the same subject appeared. Using a renormaliza- 
tion group type argument the authors 'explain' the approach to the critical 
solution and final run-away of almost critical solutions. They obtain the same 
expression for the critical exponent, which they determined numerically for 
the special case k = 1/3 of the 'Radiation Fluid'. 
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